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ABSTRACT 

Most water ice in the universe is in a form which does not occur naturally on Earth and of which only 
minimal amounts have been made in the laboratory. We have encountered this “high-density amorphous ice” 
in electron diffraction experiments of low-temperature (T < 30 K) vapor-deposited water and have subse- 
quently modeled its structure using molecular dynamics simulations. The characteristic feature of high-density 
amorphous ice is the presence of “interstitial” oxygen pair distances between 3 and 4 A. However, we find 
that the structure is best described as a collapsed lattice of the more familiar low-density amorphous form. 

These distortions are frozen in at temperatures below 38 K because, we propose, it requires the breaking of 
one hydrogen bond, on average, per molecule to relieve the strain and to restructure the lattice to that of 
low-density amorphous ice. Several features of astrophysical ice analogs studied in laboratory experiments are 
readily explained by the structural transition from high-density amorphous ice into low-density amorphous 
ice. Changes in the shape of the 3.07 pm water band, trapping efficiency of CO, CO loss, changes in the CO 
band structure, and the recombination of radicals induced by low-temperature UV photolysis all covary with 
structural changes that occur in the ice during this amorphous to amorphous transition. While the 3.07 jivw ice 
band in various astronomical environments can be modeled with spectra of simple mixtures of amorphous 
and crystalline forms, the contribution of the high-density amorphous form nearly always dominates. 

Subject headings: dust, extinction infrared: ISM: lines and bands — ISM: molecules — line: profiles — 
methods: laboratory — molecular data 


1. INTRODUCTION 

The frost on interstellar grains in dense molecular clouds is 
thought to consist mainly of an amorphous water ice with 
trapped impurities (e.g., Tielens & Allamandola 1987). The 
observed 3.07 gm absorption band is usually best compared to 
the OH stretch vibration band of impure water ice, vapor- 
deposited onto inert surfaces at temperatures less than 30 K 
(Hagen, Tielens, & Greenberg 1981 ; Schmitt et al. 1989). The 
same, or similar, ices are thought to constitute the bulk of 
comets (Mumma, We^ssman, & Stern 1993; Kouchi et al. 
1994). The 3.07 gm ice band has been studied in laboratory 
experiments of interstellar ice analogs as a function of tem- 
perature, deposition conditions, and impurity content (Hagen 
et al. 1981; Sceats & Rice 1982. Hagen et al. 1983; Tielens, 
Hagen, & Greenberg 1983; Schmitt, Ocampo, & Klinger 1987). 
In addition, considerable effort has gone into measuring the 
changing shapes of infrared absorption bands of trapped impu- 
rities, the diffusion and recombination of radicals induced by 
UV photolysis, and the outgassing properties of ice mixtures 
(Bar-Nun, Herman, & Laufer 1985; Bar-Nun et al. 1987; Grim 
& Greenberg 1987; Laufer, Kochavi, & Bar-Nun 1987; Sand- 
ford et al. 1988; Schutte 1988; Schutte, Allamandola, & 
Sanford 1993). 

The observed characteristics of vapor-deposited impure 
water ices are a function of the arrangement of the water mol- 
ecules in the solid, which can be studied directly by electron or 
X-ray diffraction. Fourier transformation of the observed dif- 
fraction patterns yields a statistical distribution of distances 
between the oxygen atoms in the water matrix (the oxygen- 
oxygen radial distribution function), which uniquely character- 
izes its structure. By applying electron diffraction to dynamic 
systems, we find that vapor-deposited water ice exists in three 
amorphous forms, apart from two crystalline states 
(Jenniskens & Blake 1994). 


This paper characterizes the two low-temperature amor- 
phous forms of vapor-deposited water ice, low-density amor- 
phous ice (I a l) and high-density amorphous ice (I a h), and the 
behavior of these ices in the temperature range 15-120 K. An 
amorphous form of water ice was first described early in this 
century. Diffraction studies were performed by Burton & 
Oliver (1935). This familiar low-density form is made when 
water vapour is deposited on an inert surface at 77 K. I a l has a 
density of 0.94 + 0.03 g cm -3 (Narten, Venkatesh, & Rice 
1976; Mayer & Pletzer 1986) and has been studied extensively 
because of its resemblance to liquid water (Olander & Rice 
1972; Sceats & Rice 1982; Pohorille et al. 1987; Sciortino, 
Geiger, & Stanley 1992), and because it can be used to replace 
liquid water as an embedding material for biological tissues in 
electron microscopy (Heide 1984; Heide & Zeitler 1985). 

High-density amorphous ice was first reported by Narten et 
al. (1976). X-ray diffraction patterns of water vapor deposited 
at 10 K on a single crystal copper substrate showed the pres- 
ence of oxygen atoms at distances between those of the first 
and second nearest neighbor at the interstitial site of the I a 1 
network, which results in a higher value for the calculated local 
density of oxygen atoms, i.e., 1.1 ±0.1 instead of 0.94 g cm 3 . 
For this reason, Narten et al. introduced the names of “high 
density” (I a h — our notation) and “low density” (I a l) amor- 
phous ice. Later experiments failed to reproduce the high- 
density form on substrates other than single-crystal copper, 
and it was thought that the high-density amorphous form was 
an artifact of epitaxial growth on an oriented single crystal 
(Sivakumar, Rice, & Sceats 1978; Sceats & Rice 1982), depend- 
ing, perhaps on a nucleation step (Sivakumar et al. 1977). 
Recently, it was shown that high-density ice also forms during 
vapor deposition on amorphous carbon films (Heide 1984; 
Jenniskens & Blake 1994). However, it has not been shown 
that the resulting ice structure is the same as that of Narten 
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et al. (1976). This high-density amorphous ice is potentially 
important in astrophysics, because infrared spectra of low- 
temperature deposits provide the best fit to the 3.07 ^xm ice 
band (Hagen et al. 1981). In the astrophysical literature, the 
high-density form has only briefly been referred to before by 
Hagen et al. 

Another form of high-density amorphous ice can be made by 
pressure melting, and is similarly referred to as high-density 
amorphous ice (hda or HDA). When cooled to 77 K, crys- 
talline hexagonal ice, I h , undergoes a transition into hda at 1-2 
GPa pressure (Mishima, Calvert, & Whalley 1984, 1985). The 
resulting ice, when recovered at 1 atm, has a density of 
1.17 ±0.03 g cm -3 (determined by weighing in liquid 
nitrogen). Theoretical estimates of the transition pressure as a 
function of temperature are given in Tse & Klein (1990) and 
Poole et al. (1992). HDA can exist metastably for some time at 
low pressure and 77 K. When this pressure-induced, high- 
density amorphous ice is warmed at a rate of 0.2-3 K 
minutes -1 , the ice transforms irreversibly into a low-density 
0.94 ± 0.02 gem -3 amorphous form (Ida or LDA) in the range 
105-123 K, releasing energy in the process (Mishima et al. 
1984, 1985; Handa, Mishima, & Whalley 1986; Bizid et al. 
1987; Tse & Klein 1987; Hemley, Chen, & Mao 1989). The 
radial distribution function of hda (Bosio, Johari, & Teixera 
1986; Bizid et al. 1987) differs significantly from that of the I a h 
of Narten et al. (1976). It has been argued that pressure- 
induced hda has more hexagonal order than vapor-deposited 
I a h (Tse 1992). 

In this paper, we confirm the existence of the vapor- 
deposited, high-density amorphous form I a h from electron 
microscopy studies (§ 2). We measured the radial distribution 
function (rdf) of I a h and the gradual change of the diffraction 
pattern during the I a h — ► I a 1 transition (§ 3, § 4). The rdf is 
compared to that of pressure-induced hda and similar data by 
Narten et al. (1976). The distinct physical properties of I a h (§ 5) 
and the conditions necessary for formation during vapor depo- 
sition (§ 6) are determined by reinterpreting experimental data 
reported in the literature. In § 7, new molecular dynamics cal- 
culations are presented that attempt to model the ultrastruc- 
ture of high-density amorphous ice. The theoretical and 
experimental results are discussed in § 8. In § 9, it is argued 
that high-density ice is ubiquitous in space, to an extent that 
this form is probably the most abundant form of ice in the 
universe. 

2. CHARACTERIZING THE STRUCTURE OF AMORPHOUS ICES 

The electron diffraction studies employ a Hitachi H-500H 
transmission electron microscope (TEM), which was adapted 
for vapor-depositing water molecules in situ onto a thin 
amorphous carbon film subtrate held at a temperature of 14 K 
(Blake & Palmer 1991). The structure of the ice film is charac- 
terized by selected area electron diffraction (SAED). Morpho- 
logical features of the ice can be viewed in either bright-field 
mode (strongly diffracting regions of the film appear dark) or 
dark-field mode (diffracting regions appear bright). 

Degassed distilled water vapor is deposited from a reservoir 
at 14 torr and room temperature into the microscope at base 
pressure of ~3 x 10 -7 torr, until the partial H 2 0 pressure 
increases to ~4 x 10 6 torr. A film of ice is grown by opening 
a shutter that covers the sample for 5-20 s (rate of growth 6-36 
jun hr -1 ), while monitoring the diffraction pattern. The 
resulting film is about 0.05 /xm thick, smooth, and homoge- 
neous. Bright-field imaging reveals no sign of the crystallinity 


reported by Olander & Rice (1972), even at high magnification. 
There is also no indication of needle- or stalactite-like surfaces 
as reported by Laufer et al. (1987). 

Diffraction patterns from the ice film are obtained from 
volumes of material of about 0.7 /xm 3 . A parallel eletron beam 
of 3 firm in diameter is directed through the sample to form the 
diffraction pattern, which is focused onto an imaging device at 
the back focal plane of the microscope. The resulting diffrac- 
tion pattern is a Fourier transform of the statistical distribu- 
tion of distances between atoms in the ice. Oxygen atoms are 
the principal electron scatterers (scattering by hydrogen atoms 
is negligible). The diffraction pattern is electronically recorded 
with an instrumental resolution of A k = 0.08 A -1 , sampled at 
0.03 A 1 pixel ! , and the patterns are stored at a rate of 1-3 
per minute during a typical experiment. Temperatures can be 
recorded with an absolute accuracy of ± 2 K and relative accu- 
racy of about 0.2 K. The warming rate in these experiments 
varied from 2.3-5.2 K minutes -1 and was typically 2.7 K 
minutes - *. 

The two-dimensional circular diffraction patterns are reduc 
ed by folding to a one-dimensional radial intensity distribution 
(Fig. 1). The x-axis, the momentum transfer coordinate 
k — (4tt/A) sin is calibrated using crystalline gold (a 0 = 
4.0786 A) deposited on a small part of the substrate as an 
internal standard, which was measured at the beginning and 
end of each experiment. The relatively small intensity contrib- 
uted by the amorphous carbon substrate is removed by a back- 
ground subtraction algorithm. 

Low electron dose imaging techniques are employed to 
minimize sample degradation during viewing. Electron beam 
damage does not change the diffraction patterns on the time- 
scale of each measurement, and different parts of the ice layer 
are exposed during each subsequent measurement. The tech- 
nique offers a powerful means of measuring the structural 
properties of water ice films in situ in vacuum. 

Fourier transformation of the diffraction patterns yields the 
rdf of distances between oxygen atom pairs (Fujime 1966; 



Fig. 1. — Representative electron diffraction patterns of high-density (I a h) 
and low-density (I a l) amorphous ice, at 14 and 77 K, respectively, after sub- 
traction of a contribution by the amorphous carbon substrate. The dashed line 
is an example of the chosen background function in the derivation of the radial 
distribution function (7 bg in eq. 1). The first and second diffraction maxima are 
labeled 44 a ** and “ b ” and are referred to in Fig. 3a-3c. 
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Moine, Pelton, & Sinclair 1988) 
k = (4rc sin &)/A , 

m = UoM - h t miPfh , 0) 

rdf = (2/ji) J* /(fc) exp ( - a/c 2 ) sin ( kr)(k/r)dk , 

where / obs is the observed diffraction pattern after subtraction 
of a substrate contribution, / bg is a smooth and continuous 
function through the inversion points of the pattern (a sum of 
three power laws; Fig. 1), /? is a scaling factor taken such that 
the intensity is — 1 for small k after subtraction of the back- 
ground (at onset of first diffraction peak), and f Q is the electron 
scattering function of oxygen. The parameter a, in the damping 
factor in the Fourier transform, is taken such that exp 
( — a/c^ ax ) = 0.1. The latter reflects the rapid decrease of the 
electron beam intensity away from the center of the beam and 
the limited range of k that is measured. 

Due to limitations in the dynamic range of our electronic 
imaging system, we cannot derive an rdf for every recorded 
pattern. In order to follow the changes in the structural state, 
we measure the position, height, and width of each peak in the 
diffraction pattern by fitting Gaussian profiles to individual 
maxima. The Gaussian profiles usually produce a good fit to 
the peaks, which is the only and sufficient reason for choosing 
this particular functional form. 

3. THE RADIAL DISTRIBUTION FUNCTION OF I a h AND H 

When water ice is formed by vapor deposition at 14 K, it 
yields a distinct diffraction pattern when compared to that of 
the more familiar low-density amorphous ice at 77 K (Fig. 1). 
Most notably, the first diffraction peak is shifted to larger 
values of k. 

The resulting radial distribution function of the low- 
temperature deposit (Fig. 2), which Narten et al. (1976) calls 
“high-density amorphous ice” (I a h — our notation), differs 
principally from the higher temperature deposit, I a l, in having 
an additional maximum between 3 and 4 A. 

It is illustrative to compare the radial distribution functions 
of I a h and I a l to those of vapor-deposited amorphous ice at 10 
K (Narten et al. 1976), and those of pressure amorphosed hex- 
agonal ice (Bizid et al. 1987). Figure 2 shows that there is a 
striking similarity between our rdf of vapor-deposited I a h and 
that of the pressure-induced D 2 0 ice by Bizid et al. (1987). The 
high-density amorphous form is characterized, in both cases, 
by a splitting of the second nearest neighbor peak, such that 
one component lies between the first and second nearest neigh- 
bor peaks of the low-density amorphous form. The diffraction 
maxima in I a h are less pronounced, suggesting that I a h is less 
well ordered than BizicTs hda. This apparent reduction in ord- 
ering may be due to the thinness of the vapor-deposited ice, 
differences in annealing times between the two ices, exposure to 
high fluxes of X-rays, or more microporosity in the vapor- 
deposited ice (that is, defects in the random hydrogen-bonded 
network with typical cavity sizes of radius 1 A < r < 10 A). 

The agreement with the pressure-induced hda is better than 
with the vapor deposit of Narten et al. (1976). The first nearest 
neighbor distance of both I a h and I a l is 2.8 1(T 0.05/ — 0.01) A, 
in good agreement with the hda of Bizid et al. (1987) (2.80 A) 
and Bosio et al. (1986) ( ^2.82 A), while Narten et al. (1976) 
report a slightly smaller value at 2.76 A. The second nearest 
neighbor in our I a h is broad and is split into two peaks at 


3.70 A and 4.70 A, similar to the hda of Bizid et al., 3.70 and 
4.64 A; and that of Bosio et al., 3.65 and 4.45 A, respectively. 
However, neither our, nor Bizid’s, radial distribution functions 
show a sharp spike at the “interstitial” position as in the 
Narten et al. (1976) data. 



Fig. 2. — Fourier transform of the diffraction pattern. These oxygen-oxygen 
distance radial distribution functions (rdf's) depict the probability of finding 
another oxygen atom at given radial distance. The rdf of low-density ice is 
displaced by —1.0 for clarity. From top to bottom: (a) rdf derived from our 
electron diffraction data of vapor-deposited ice at 15 K. (top) and after 
warm-up to 77 K. ( bottom ); ( b ) X-ray diffraction data of the pressure-induced 
high- and low-density ice at 77 K. by Bizid et al. (1987) directly after pressure 
relief and after a brief warming to 140 K and subsequent cooling, respectively; 
and (c) as derived by Narten et al. (1976) from X-ray data of vapor-deposited 
water at 10 K (I, h) and of a vapor deposit at 77 K (I a 1). 
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also be identified indirectly. We have reinterpreted experimen- 
tal data reported in the literature, in light of the previous 
results, and suggest that the high-density amorphous form can 
be identified by the shape of the 3.07 /xm band (Hagen et al. 
1981; Mayer & Pletzer 1985), the shape of the dangling OH 
bond absorption bands (Venkatesh, Rice, & Bates 1975; Siva- 
kumar et al. 1978), the presence of CO impurities (Sandford et 
al. 1988; Schmitt, Greenberg, & Grim 1989), outgassing behav- 
ior (Laufer et al. 1987; Kouchi 1990), and the diffusion of rad- 
icals (Schutte 1988). 

Figure 4 summarizes the deposition conditions for the 
experiments cited above. It is concluded that I a h forms readily 
when water is condensed slowly (<100 /xm hr 1 ) at low 
enough temperature ( T <40 K) on any surface with good 
thermal conductivity (single-crystal copper, amorphous- 
carbon film, aluminum, gold-coated copper, Csl window). 
Vapor deposits at 40-70 K form mostly I a 1 (see, e.g., Sandford 
etal. 1988). 

There is a maximum and a minimum condensation rate at 
which I a h can form. Deposition rates close to the maximum of 
100 /xm hr -1 do result in a structure with more I a l character. 
For example, Palumbo & Strazzulla (1993) measured a notice- 
able change in the position of the C—O stretch in an impure 
water ice for deposition rates above 1 /xm hr' 1 . For a given 
temperature, there is, potentially, a minimum deposition rate, 
below which the water molecules have time to reorient them- 
selves before being incorporated in the I a h lattice by sub- 
sequent condensing molecules. Kouchi et al. (1994) measured 
this effect for the formation of the cubic crystal form in the 
low-density amorphous domain at very low deposition rates, 
and found a temperature dependence, which, to us, suggests a 
characteristic A H — 44 + 7 kJ mol ” 1 activation enthalpy. This 
value is the same as that for the crystallization of I a 1 during 
warm-up (Jenniskens & Blake 1994). An extrapolation of the 
measured critical flux at 96-1 10 K (Kouchi et al. 1994) down to 



deposition rate (micron/hr) 

Fig. 4. — Deposition conditions that are favorable for the formation of I a h 
(diamond) and I 4 1 (cross), as inferred from five deposition experiments at 14 K 
and from direct and indirect indicators in 23 literature studies. A combination 
of diamond and cross indicates a weak high-density amorphous ice signature 
only. The one at T = 15 K is a thick 1 mm ice layer by Venkatesh, Rice, & 
Bates (1975). The dashed lines are approximate boundaries discussed in the 
main text. 


lower temperature suggests that, at realistically low deposition 
rates, no cubic crystalline ice will be formed in the high-density 
regime below 40 K. However, a lower activation energy may 
apply for the reorientation of the water into an open I a I struc- 
ture. The horizontal line in Figure 4 shows this deposition rate 
dependence, if the characteristic activation enthalpy is 12 kJ 
mol ' 1 (see § 8.2). The result suggests that this effect plays no 
significant role in the laboratory experiments, but may be rele- 
vant in case of extremely low deposition rates in some astro- 
physical environments. 

The film thickness has an effect as well. Most experiments 
concern layers of 0.01-10 /xm in thickness. Low-temperature 
vapor deposits as thick as 1000 /xm sometimes do (Narten et al. 
1976), and sometimes do not, show I a h signature (Venkatesh et 
al. 1975). Perhaps this is because in thicker films the heat of 
fusion from impinging molecules is not conducted fast enough 
to prevent local restructuring. There is also some indication 
that for very thin layers (<0.002 /xm), substrate surface struc- 
ture can affect the structural growth during deposition (Moore 
et al. 1994) and, therefore, the presence of a surface film on the 
substrate, and the existence of single-crystal orientations of the 
substrate, may affect the structure of thin deposits through 
epitaxy or nucleation (Sceats & Rice 1982). 

Finally, we note that the I a h->I a l transition does not 
always unfold the same way in vapor-deposited layers. For 
example, Hagen et al. (1981) studied a 0.5 /xm thick water ice 
layer, vapor-deposited at 10 K. We assume that the observed 
shift of the peak frequency and width of the infrared O—H 
stretching band is due to the I a h — ► I a I transition, and find that 
the transition had an onset shortly after the beginning of the 
warm-up at about 20 K, which is at lower temperature than in 
our data. The transition was complete at about 70 K, as in our 
data. On the other hand, the 1000 /xm thick ice layer of Narten 
et al. (1976) still had a recognizable “interstitial” peak in the 
rdf after warming to 77 K. 

7. MODEL CALCULATIONS 

Molecular dynamics (MD) calculations of water structure, 
that is, solutions of the equations of motion for an ensemble of 
point masses with the intermolecular potentials of water, go 
back two decades (for reviews see Angell 1983 and Dore & 
Teixeira 1991). The resulting structures are best described as an 
imperfect random network of hydrogen bonds. The models 
successfully describe many features of the measured structure 
and physical properties of liquid water and low-density ice I a 1. 

7.1. Vapor Deposition 

MD simulations of water vapor deposition at 10 K on small 
water clusters (Buch 1992) or on films (Wilson et al. 1995) do 
not yield I a h. Instead, the resulting network is found to contain 
a substantial fraction of molecules with hydrogen bond coordi- 
nation different from 4 (1, 2, 3 and 5) and an abundance of 
dangling OH bonds. This results in a porous structure with 
many microvoids and a very large effective surface area. Buch 
(1992) ascribed this behavior to the preferred attachment of 
new molecules to dangling atoms on the surface (unbonded H 
atoms) and to O atoms with only one hydrogen bond, thus 
propagating the formation of surface protrusions. MD calcu- 
lations are limited by the amount of relaxation that can be 
allowed after a molecule has collided with the surface. The 
ensemble can usually only be followed for tens of picoseconds 
before allowing the next molecule to hit the surface. In effect, 
this corresponds to deposition rates of 10 8 -10 9 /xm hr' 1 . The 
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Fig. 5 —Model of water ice. A 10 x 10 x 10 A 3 sample from a MD simula- 
tion of (a) high-density amorphous ice ( b ) low density amorphous ice, and {c) 
cubic crystalline ice. These stick diagrams are randomly chosen from one final 
configuration of a 340 molecule model and show the position of the oxygen 
atoms only. The van der Waals radius of the water molecules is larger than the 
size of the symbols. Solid lines connect the oxygen atoms that are less than 3.0 
A apart. For comparison purposes, the cubic ice is shown at an unfavorable 
angle with respect to the crystal lattice planes and axis of symmetry. 


molecular velocities within the cluster are rescaled, after each 
sticking trajectory, by a constant factor to represent cooling 
back to the initial temperature of 10 K. However, this quen- 
ching may not result in the same final state as a fully relaxed 
network. In addition, the water molecules are considered rigid 
and certain energy transfer pathways are not considered in the 
calculations, including quantum mechanical tunneling and 
energy transfer to intramolecular O — H stretch vibrations. 
Note that the energy of binding is of similar magnitude as 



Ial 


Fig. 5b 

the energy of an intramolecular O— H stretch vibration or 
bending mode. 

7.2. Pressure Melting 

MD simulations have modeled the transition of hexagonal 
crystalline ice (I h ) into high-density amorphous ice by pressure 
(Tse & Klein 1987; Tse 1992). The process was found to be 
driven by a collapse of the open hexagonal ice I h structure. In 
the resulting oxygen-oxygen rdf, the second nearest neighbor 
peak, at 4.5 A, broadens significantly, and a new band, at 3.2 A, 
starts to emerge, reflecting a migration of water molecules into 
the interstitial sites. A number of hydrogen bonds are broken, 
creating 3-coordinated and 5-coordinated molecules. Some 
memory of I h remains: the network appears to keep a recog- 
nizable layering (stacking order) and 6-fold symmetry when 
viewed along the corresponding crystal axis (Tse 1992). 

7.3. Quenching Configurations of Liquid Water 

We have attempted to model high-density amorphous ice by 
quenching configurations, drawn from MD simulations of 
liquid water in a confined space at 300 K, instead of a crys- 
talline I h configuration, therefore starting from an amorphous 
configuration. The system consisted of 343 water molecules 
placed in a cubic simulation box, assuming periodic boundary 
conditions and imposing a density of 1.17 g cm -3 . Instead of a 
steepest descent minimization of the high-density water con- 
figurations, as in our previous model (Wilson et al. 1995), a 
more gentle quenching was attempted. Each stage of quench- 
ing consisted of a 20 ps trajectory during which the tem- 
perature was reduced by 50 K, followed by 20 ps of equili- 
bration. Most of the network changes occurred in the higher 
temperature regime above 150 K. Subsequently, at 77 K, the 
simulation box was expanded in one dimension, thereby cre- 
ating a slab of ice with two free surfaces. The structure did not 
change on a tens of picoseconds timescale when allowed to 
expand in two dimensions. 

To illustrate the disordered nature of these amorphous ices, 
Figures 5a-5c show a small 10 x 10 x 10 A cube of the high- 
density amorphous, low-density amorphous, and crystalline I c 
configurations. The planes of molecules, and the cubic and 
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hexagonal symmetries that are typical of the crystalline poly- 
morphs, cannot be recognized in the amorphous form when 
viewed from different angles. There are many 4-, 5-, and 7- 
ringed structures which create a wide range of cavities or 
“ interstitial ” sites. 

The calculated oxygen-oxygen rdf of the high-density 
amorphous ice, at 77 K, shows a broad second nearest neigh- 



Fig. 6. — Radial distribution functions, averaged from 50 such configu- 
rations of high- and low-density ice as shown in Fig. 5. The oxygen-oxygen rdf 
can be compared to those in Fig. 2: I a h has a broad second nearest neighbor 
peak, with a shoulder at 3.7 A that is not found in I a 1. The hydrogen-hydrogen 
rdf (only inte molecular distances) shows the same pattern of bond distances. 
There is no clear difference between high- and lowndensity ice in the oxygen- 
hydrogen rdf (again only intermolecular distances). 


bor oxygen peak, with some indication of splitting into two 
components (Fig. 6). The peak is not as pronounced as in the 
collapsed I h model of Tse & Klein (1987), but is clearly different 
from the rdf of the low-density ice model. Also, this peak spans 
the region between the “ interstitial ” and second nearest neigh- 
bor peak in the observed rdf of I a h (Fig. 2). 

These models of I a h and I a l were investigated further by 
calculating the distribution of coordination numbers and the 
distribution of cavity sites in the ice, as in Pratt & Pohorille 
(1992). The coordination number (Table 2) is defined in two 
ways: either by the number of oxygen atoms within a sphere of 
3.2 A, or by the presence of a hydrogen bond stronger than — 8 
kJ mol -1 . The number distribution of hydrogen bonds shows 
that (1) there are no interstitial water molecules with zero 
hydrogen bonds; (2) the most preferred coordination is 4; (3) 
the distribution of low coordination numbers less than 4 does 
not differ between the high and low-density amorphous ice 
models; (4) in high-density amorphous ice, there is a large 
percentage of molecules with five, six, and even seven oxygen 
atoms located within 3.2 A, that are not found in low-density 
ice; (5) the increased percentage of short oxygen-oxygen dis- 
tances does not correspond to a large number of 5-coordinated 
hydrogen bonds, although the percentage of 5-coordinated 
water is somewhat higher in high-density ice. 

Figure 7 shows the probability of placing a rigid sphere of 
given radius at a random location in the ice structure (allowing 
for the Van der Waals radius of the water molecules). Large 
open spaces are conspicuously absent in high-density amor- 
phous ice, relative to low-density amorphous ice. For radii 
larger than 0.3 A, there is a fast decline of the relative probabil- 
ity of placing a cavity of given radius at given position. In the 
model of high-density amorphous ice, cavities larger or equal 
to the largest cavity in the cubic crystalline lattice are highly 
improbable. 

8. DISCUSSION 

8.1. A Model of High- Density Ice 

Two models have been proposed for vapor-deposited high- 
density ice: (1) unbonded water molecules are trapped in 
interstitial sites in slightly distorted ice I cages, in which 45% of 
the cavities are filled with interstitial water, or 19% of all water 
molecules are in cavity positions (Narten et al. 1976); and (2) 
high-density amorphous ice is a randomized version of ice II or 
III, high-pressure crystalline polymorphs of water ice which 

TABLE 2 

Percentile Distribution of Oxygen Coordination Numbers* 

0-0 < 3.2 A H Bond < -SUmoE 1 

Coordination 


Number hda Ida I c hda Ida I c 


0 0.00 0.00 0.00 0.26 0.18 0.00 

1 0.00 0.00 0.00 4.9 4.4 0.00 

2 0.02 0.13 0.00 17.1 16.3 0.00 

3 2.5 6.9 0.00 28.3 28.4 0.00 

4 58.8 86.4 100.00 43.9 48.7 100.00 

5 30.3 6.4 0.00 5.4 2.00 0.00 

6 7.3 0.15 0.00 0.08 0.00 0.00 

7 0.06 0.003 0.00 0.00 0.00 0.00 


* Based on the number of oxygen atoms within a radius of 3.2 A and the 
number of those oxygen atoms that have a hydrogen bond strength larger 
than - 8 kJ moP 1 (statistical error is ± 2 in the least significant digit). 
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Fig. 7. — The probability of placing a rigid sphere of given radius at a 
random location in the free space of the water ice model, allowing for the van 
der Waals radius of the water molecules. The three curves are for the models of 
the two low-temperature amorphous forms and the crystalline cubic lattice. 
The high-density amorphous form has a significant lack of the large cavities 
found in the low-density amorphous form. 


have a similar mean density (Narten et al. 1976; Sceats & Rice 
1982). 

These models are highly unlikely in a distorted network. 
Some bonding of the interstitial water to other water molecules 
is expected. The presence of a distorted form of one of the 
high-pressure crystalline polymorphs of ice is also, in our 
opinion, not reasonable, because either these polymorphs need 
a form of proton ordering (preferred orientation of the water 
molecules), or they are outside their stability regime. 

Following Tse & Klein’s interpretation of the difference 
between hda and Ida (Tse & Klein 1987; Tse 1992), we propose 
an alternative model: the high-density amorphous form is a 
collapsed low-density lattice, characterized by the presence of 
oxygen-oxygen distances between 3 and 4 A. Arguments for 
using the term “collapsed” include that the number density of 
the large cavity sites is significantly reduced in high-density 
amorphous ice compared to the low-density ice. Also, there are 
more cases of 5 oxygen atoms within 3,2 A of another oxygen 
atom, without a similar increase in the number of 5- 
coordinated hydrogen bonds. 

This raises the question of how a collapsed lattice is formed 
during vapor deposition. Indeed, MD model calculations typi- 
cally result in the formation of an open I a l lattice. There is a 
possibility that ice is deposited in an I a h form because nature 
offers a faster redistribution of the energy of binding (some 30 
kJ moP 1 on average) than the molecular dynamics calcu- 
lations allow for (e.g., because in the MD model the water 
molecules are rigid). Alternatively, the ice may deposit itself in 
an open I a l form, but a collapse to I a h occurs below the 
binding site upon release of binding energy, analogous to the 
injection of energy in the ice network, as in charged particle 
bombardment and photolysis (§ 4), while the MD calculations, 
again, do not allow a fast enough dissipation of energy for the 
distortion to be frozen in. 

There is no single structure which characterizes I a h. Instead, 
a wide range of configurations is possible, depending on the 
number of collapsed cages and the amount of relaxation 
toward an I a 1 structure. Such relaxation apparently causes the 


less distorted cages to be restructured more easily than the 
more heavily distorted structures, with water molecules at an 
“interstitial” position in an ice I a l lattice. The difference 
between hda and I a h, if any, is that hda may have a more 
recognizable hexagonal stacking order than I a h. 

8.2. I a h^> I a l and l a l-> I a h 

The presence of a well-defined transition between one 
amorphous form and the other suggests that there is an activa- 
tion energy barrier that prevents the transition at temperatures 
below 38 K in a vapor deposit, and significantly delays the 
transition below 70-80 K in electron-irradiated I a l or I c . The 
nature of this barrier determines the physical difference 
between the two amorphous forms. 

We speculate that the I a h — ► l a 1 transformation requires the 
breaking of one hydrogen bond, on average, per molecule in 
the I a h structure. An argument for this is that the temperature 
range of 70-90 K that discriminates the two regimes is at half 
the temperature needed to transform amorphous ice into cubic 
ice. The latter transition occurs at 143-155 K and has a char- 
acteristic activation enthalpy of A H a = 44 + 2 kJ mol 1 
(Jenniskens & Blake 1994), which is equivalent to breaking 2 
hydrogen bonds on average per molecule. The onset of the 
transition will scale with the activation energy, if the kinetic 
mechanism that drives the transformation is the same. There- 
fore, the breaking of 1 hydrogen bond on average per molecule, 
with an activation energy of Mi a = 22 kJ moK \ corresponds 
to a temperature of about 72-78 K, in agreement with the 
boundary of the metastability zones of high and low-density 
amorphous ice (70-90 K). 

The reordering from I a h to I a l starts at 38 K, in a gradual 
warm-up, and proceeds gradually. The onset at 38 K suggests 
that the activation energy is lowered by some mechanism, and 
the gradual transition indicates that the effect is of a range of 
activation energies. This probably reflects the distribution of 
binding states in the ice. Indeed, MD models of vapor- 
deposited low-density ice (Buch 1990; Zhang & Buch 1990; 
Buch 1992) show many defects in this network, causing a range 
in bond coordination. The model low-density ice has a mean 
potential energy per molecule of 20.9 kJ mol \ averaged over 
time and all molecules in a cluster, but a broad distribution 
extends from about 7 to 28 kJ mol 1 and peaks at 17.5 and 
23.5 kJ mol 1 . Alternatively, it has been shown that 5- 
coordination can lower the activation enthalpy by about 10 kJ 
mole' 1 (Sciortino et al. 1992). High-density amorphous ice is 
particularly rich in 5-coordinated water. 

Exposure to the electron beam at temperatures below 70 K 
causes the ice to flow, although the I a h structure is preserved. 
The low viscosity of high-density ice suggests that it can easily 
restructure its lattice, like liquid water. In liquid water, 5- 
coordination has been invoked to account for the easy 
restructuring, because it allows a lowering of the activation 
enthalpy (Sciortino et al. 1992). Five-coordination is less 
important in low-density ice, and absent in cubic ice, in agree- 
ment with the observed higher viscosity. 

The mechanism of electron beam processing of I c into I a h is 
not yet well understood. Elastic collisions of electrons in the 
bulk are not important for displacements (Heide 1984) and, 
hence, for the formation of I a h. At a given dose, only about 
2 x 10' 4 of all molecules are displaced. Ionization processes, 
in which a water molecule is dissociated (H 2 0->H + OH), 
can be significant. For the given radiation dose, about 30% of 
all water molecules are dissociated. However, other processes 
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may dominate (for example, the electronic or vibrational exci- 
tation of the inter- and intramolecular hydrogen bonds). 
Studies of the influence of photons on this process may reveal 
new information. Note, for example, that photons that excite a 
single O-H-O bending or a single O-H stretching vibration 
(6.27 and 3.07 p m) carry as much energy as a typical hydrogen 
bond (19 and 39 kJ mol -1 , respectively). Visible and near- 
infrared photons carry more energy than this and may contrib- 
ute to the structural transformation if absorbed. 

8.3. Time-Temperature Diagram 

It is reasonable to assume that activation enthalpy is a local 
phenomenon and activation enthalpies at other sites in the 
lattice are not much affected by a local restructuring. Given a 
range of activation enthalpies that determine the onset and the 
completion of the high- to low-density transition, a time- 
temperature diagram can be constructed (Fig. 8), according to 
which the induction time (At) of each phase of the transition 
behaves according to (Schmitt et al. 1989) 

At(s) = v 0 exp (AHJRT) . (2) 

The kinetics of all transitions in water are governed by the 
intermolecular vibration (v 0 = 9x 10 -15 s) or bending (v 0 = 
2.1 x 10“ 14 s) timescales of water. The uncertainty in the time 
axis of Figure 8 (the vertical shift) is about a factor of 4, 
depending on the definition of the onset. 

If AH a is of order 22 kJ mol -1 , it follows that high-density 
amorphous ice can persist on a timescale of 10 5 yr at tem- 
peratures below 22 K, even in the absence of irradiation or 
charged particle bombardment. Similarly, the high-density ice 
can persist below temperatures of 19 K on the solar system 
timescale of 4.5 x 10 9 yr. 



Fig. 8. — Time-temperature diagram, based on conjecture, showing the time 
it takes for the ice to start and complete the transformation from I, h to 1, 1 
after being warmed to a given temperature. The assumptions made are 
explained in the main text. Markers at 4.5 x IO 9 and 10 3 yr indicate the typical 
timescale of water-ice evolution in the solar system and in the dense interstellar 
medium, respectively. 


9. ASTROPHYSICAL IMPLICATIONS 
9.1. The 3.07 pm Ice Band on Interstellar Dust Grains 

Ice in the interstellar medium is usually characterized by the 
strength and shape of the 3.07 pm absorption band in the 
infrared (Knacke, Cudaback, & Gaustad 1969; Gillett & 
Forrest 1973; Merrill, Russell, & Soifer 1976; Hagen et al. 
1981; Whittet et al. 1983; Bult, Greenberg, & Whittet 1985; 
Baratta & Stazzulla 1990; Smith, Sellgren, & Brooke 1993). 
This band is due to the O— H stretch vibration of water. A 
unique assignment of shoulder and peaks of bands to vibra- 
tions is not possible because intra- and intermolecular coup- 
ling of O—H oscillators are nearly equal in strength (Haas & 
Hornig 1960; Bergren et al. 1978; McGraw et al. 1978). Inco- 
herent inelastic neutron scattering gives a density of vibra- 
tional states that shows a strong difference between 
pressure-induced high- and low-density amorphous ices (hda 
and Ida) (Klug et al. 1991 ; Tse 1992). The strongest differences 
are in the longitudinal and transverse optic branches of the 
dispersion relation. The leading edge of the librational band 
begins at lower frequency in hda ice and rises more slowly. 
These differences are manifested in the band shapes of infrared 
absorption bands. 

There is a clear difference in the 3.07 ^m band shape associ- 
ated with vapor-deposited high- and low-density amorphous 
ice, I a h, and I a l. If the laboratory spectra of low-temperature 
deposits by Hagen et al. (1981) are indeed those of I ft h, as we 
believe they are, it follows that the I a h band contains addi- 
tional absorption at high vibrational frequency. When the ice 
is warmed at a rate of less than 5 K minutes” 1 , the peak 
wavelength changes from 3.04 to 3.08 pm , and the feature 
sharpens from 310 to 230 cm -1 in the temperature range 
20-70 K. The 45 pm band is a less sensitive probe of the 
I a h -► I a 1 transition (Klug et al. 1991 ; Moore & Hudson 1992). 

Smith, Sellgren, & Tokunaga (1989) found that some astro- 
physical objects show ice-band profiles that are best fitted by a 
combination of ice deposited at low temperature and warmed 
to 23, 77, and 150 K. We now interpret these deposits as high- 
density amorphous ice I a h, low-density amorphous ice I a 1, and 
cubic crystalline ice I c , respectively. For example, spectra from 
the Becklin-Neugebauer object are best fitted with spectral 
ratios of I a h:I a l:I c equal to 3.1 : 1.2: 1, while Mon/IRS 3 is best 
fitted by ratios of 5.3 : 2.5 : 1 and GL 259 1 by ratios of 8.6 : 1 .2 : 1 . 
The 23 K ice band provides a satisfactory fit to the overall 
shape of most other known spectra. The conclusion must be 
that ice frozen on grains in the interstellar medium is prin- 
cipally in the high-density amorphous form. 

The frost on interstellar grains constitutes the largest 
reservoir of ice in the universe — in our Galaxy alone about 
10 37 kg (Jenniskens et al. 1993). Little is known about extra- 
solar planetary systems, hence the amount of ice in icy planets 
and comets cannot be estimated with any certainty. However, 
if all stars have an amount of icy material in the form of icy 
planets, satellites and Kuiper Belt comets similar to that in our 
planetary system, then the total adds to perhaps 10 3O -10 31 kg. 
We cannot think of other sources of ice. Hence, we conclude 
that the frost on interstellar grains is the largest reservoir of ice, 
and most ice in the universe is, therefore, in the high-density 
amorphous form. 

9.2. Formation of Icy Mantles 

The physical conditions present in molecular clouds are con- 
sistent with the presence of I a h instead of I a 1. The deposition 
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rate of water vapor in dense clouds on 10 K cold silicate grains 
can be as low as 2 x 10 -11 p m hr -1 , assuming that a typical 
ice layer of 0.02 pm is grown in 10 5 yr (Jenniskens et al. 1993). 
Extrapolation of experimental data shows that, in spite of this 
very low deposition rate, no I c crystals form at substrate tem- 
peratures below 66 K, given the activation energy involved in 
the I a 1 -* I c transition. If, similarly, some activation energy 
applies for the I a h -*> I a 1 transition, it can be calculated 
whether I a h or I a l forms. For an activation enthalpy of 12 kJ 
mol -1 , I a h, and not I a l, is formed for any substrate tem- 
perature below 30 K. Hence, evaporation and recondensation 
of icy mantles in dense cloud environments results in I a h ice 
mantles, as observed. 

Initially, however, the ice mantle is not formed by direct 
vapor deposition but by surface reactions of atomic oxygen 
and hydrogen: O + H — ► OH and OH + H -► H 2 0. In some 
environments, this initial ice layer may be preserved. One 
wonders if that ice would be structurally different from vapor- 
deposited amorphous ice. The energy released in these surface 
reactions is on the order of 1 eV, which will locally heat the 
surface. Kouchi et al. (1994) argue that the resulting tem- 
perature spike decreases rapidly with increasing distance from 
the formation point, sufficiently fast to disallow the water mol- 
ecules from setting into a crystalline structure. But does this 
prevent the formation of I a l instead of I a h? It probably does. 
We point to the analogy with electron bombardment, where 
the input of energy in the ice matrix at T < 70 K results in 
high-density amorphous ice. Hence, if the release of the binding 
energy of O— *H chemical bonds has an effect similar to that of 
the deposition of energy by electrons or protons, then the ice is 
also expected to be distorted into the I a h form during initial 
formation of the water mantle at low temperature. 

9.3. CO as a Tracer of Structural Changes in the Water Matrix 

Impurities in the water ice matrix can be used as probes of 
the structure of the water mantle (Grim & Greenberg 1987). 
The band positions of CO and C0 2 have been studied in the 
laboratory (Sandford et al. 1988; Schmitt et al. 1989; Sandford 
& Allamandola 1990; Palumbo & Strazzulla 1993). Both 
bands are readily observable in astrophysical ices. Notably the 
4.7 pm band of C=0 has been studied using ground-based 
infrared telescopes (Whittet et al. 1983; Lacy et al. 1984; Kerr, 
Adamson, & Whittet 1993; Chiar et al. 1994). 

Two absorption maxima were identified in laboratory 
studies of C0:H 2 0 1 : 20 mixtures at 10 K (Sandford & Alla- 
mandola 1988). The main absorption is at 2137 cm' 1 , which is 
shifted to lower frequency (i.e., a repulsive bond) with respect 
to the position of pure CO at 2143 cm " l . In addition, there is a 
shoulder at 2152 cm -1 (an attractive bond) in fresh deposits at 
T < 40 K, which accounts for 22%-26% of the total band 
intensity below 35 K. The latter band progressively disappears 
(irreversibly) when the ice is warmed above 40 K. The CO loss 
stops immediately when warm-up is stopped (Sandford & Alla- 
mandola 1988), which implies that the loss of CO is controlled 
by changes in the water structure and that these structural 
changes are gradual. Indeed, the temperature zone where this 
band disappears (40 60 K at 2 K minutes -1 warm-up) coin- 
cides with our measured transition I a h — ► I a l. 

The peak at 2152 cm 1 was assigned to interstitial CO by 
Sandford & Allamandola (1988), and the peak at 2137 cm 1 
was thought to be due to CO that substitutes for a water 
molecule in the water-ice network. It is unlikely that CO can 
substitute for water in a water matrix, particularly at tem- 


peratures above 60 K, and that upon warming from low tem- 
peratures, interstitial CO evolves into a substitutional site. 
Schmitt et al. (1989) and Grim (1988) assigned the two peaks to 
two types of water-CO interactions within the incompletely 
hydrogen-bonded network. While the 2137 cm 1 peak was 
thought to be due to unbonded CO, the 2152 cm 1 peak was 
assigned to C0-H 2 0 hydrogen-bonded complexes in a water 
ice. This complex has a peak at 2149 cm 1 in an argon matrix 
(Dubost 1976) and 2146.5 cm 1 in a CO ice (Hagen & Tielens 
1981). During warm-up from 10 to 40 K, the hydrogen-bonded 
complexes were thought to be preferentially replaced by 
hydrogen bonds of the normal water network. We believe this 
interpretation to be essentially correct. In light of our obser- 
vation that the 2 1 52 cm - 1 peak is linked to I a h ice, we propose 
that this peak is due to CO-H 2 0 complexes involving dangling 
OH bonds, while the 2137 cm “ 1 peak is due to CO trapped in 
“ interstitial” sites in the ice cages, noting that such cages in 
I a h have a broader range of shapes and are typically smaller 
than those in I a l. This can account for the small shift to lower 
frequency and the slight narrowing of the C=0 stretch vibra- 
tion during the I a h — ► I a I transition (Sandford & Allamandola 
1988). 

The CO peak at 2152 cm - 1 is absent in interstellar ices. This 
implies that the CO-containing ice mantles were warmed 
above 60 K on timescales less than that for water mantle 
evaporation. Indeed, most of the CO is found in more volatile 
apolar mantles. However, the water ice itself is not found in the 
I a l form. This observation suggests that I a l is transformed 
back into I a h on short timescales. The most likely transform- 
ation mechanism is photolysis. The photon flux in the dense 
cloud medium is high enough to allow a restructuring into the 
I a h form on a timescale of less than 10 5 yr, when 1 %- 10% of 
water molecules absorb a 10 eV photon. Many more lower 
energy photons are absorbed. This cycling of I a h — ► I a 1 ice and 
back can play an important role in the generation of organic 
compounds of large molecular weight through UV photolysis, 
as long as traces of CO or C0 2 remain in the ice (Jenniskens et 
al. 1993; Jenniskens & Blake 1994). The photolysis will also 
cause the ice to be less viscous and form a smooth mantle 
around the ice grains. Hence, the effective surface area for H 2 
formation is probably not much different from that of a homo- 
geneous surface. 

10. CONCLUSIONS 

The observed changes of the diffraction pattern of vapor- 
deposited water ice as a function of temperature, the presence 
of two stability regimes for electron bombardment, and the 
similar rdf compared to the pressure-induced high- and low- 
density amorphous forms, among other arguments, leave no 
doubt as to the existence of high-density amorphous ice in 
low-temperature vapor deposits. It is a naturally occurring 
form which arises when water ice is slowly (< 100 pm hr -1 ) 
vapor-deposited on any substrate of good thermal conductivi- 
ty at temperatures T < 30 K. 

The radial distribution function of high-density amorphous 
ice, obtained by Fourier transformation of the diffraction pat- 
terns, has oxygen-oxygen distances between 3 and 4 A, between 
the first and second nearest neighbor distance in low-density 
amorphous ice. The shape of the rdf is similar to that of 
pressure-induced hda reported by Bizid et al. (1987). 

MD calculations can reproduce this 3-4 A feature, if the 
water ensemble is allowed to relax in a confined space. The 
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feature is not due to unbonded interstitial water molecules. 
Four-coordination is still preferred, although in comparison to 
the low-density amorphous ice model, there is a somewhat 
larger number of five-coordinated water molecules. The most 
striking feature of the high-density amorphous structure is the 
small probability of finding cavities as large as those found in 
the two crystalline forms of ice, while low-density amorphous 
ice has a cavity size distribution more characteristic of an open 
crystalline structure (Fig. 7). From these results, we conclude 
that the high-density amorphous form is best described as a 
(partially) collapsed low-density ice structure. 

There is an irreversible first-order transition from high- to 
low-density amorphous ice. Below 70-80 K, there is a range of 
activation energy barriers that hinder restructuring. Several 
processes can transform the low-density amorphous form (or 
the crystalline cubic polymorph) back into high-density 
amorphous ice in this regime. These include electron and 
proton bombardment, pressure melting, and possibly UV pho- 
tolysis. From similarities between the resulting radial distribu- 
tion functions, and from results of MD simulations of 
pressure-induced amorphous ice in the literature, we infer that 
all of these processes cause some form of collapse of the low- 
density amorphous ice structure. During vapor deposition, the 
collapse of the water lattice may be a consequence of a rapid 
dissipation of the energy of binding before impinging mol- 
ecules can reorient to form low-density amorphous ice. Alter- 
natively, the ice may form a low-density amorphous ice 
structure initially, but then be distorted by the release of 
binding energy (analogous to distortions induced by electron 
irradiation). 

The nature of the energy barrier determines the physical 
difference between the two amorphous forms. We propose that 
the restructuring of high-density ice into low-density ice 
requires the breaking of one hydrogen bond, on average, per 
molecule. The wide range of activation energies present at dif- 
ferent locations in the vapor-deposited high-density amor- 
phous ice may be responsible for the early onset. Alternatively, 
high-density ice does have a relatively large fraction of five- 
coordinated water molecules, and there has been a suggestion 
that five-coordination is responsible for the lowering of the 
activation enthalpy needed for the rapid restructuring in liquid 
water. 

Given the astrophysical conditions in the dense interstellar 
medium, we argue that vapor deposition on interstellar grains 


results in a water frost with high-density amorphous structure. 
We surmise that water mantles formed through the recombi- 
nation of H and OH radicals result also in a high-density ice 
structure. A brief warming above 40 K, however, will lead to a 
low-density amorphous form in the absence of sufficient irra- 
diation processes. The structural transition from one amor- 
phous form into another can play an important role in 
dynamic processes like radical diffusion and outgassing. 

Ground-based observations of the O— H stretch vibration 
at 3.07 pm in a variety of astrophysical environments have 
been interpreted in terms of a composite of ice deposited at low 
temperature and warmed to 23, 77, and 150 K. We now inter- 
pret these deposits as high-density amorphous ice (I a h), low- 
density amorphous ice (I a l), and cubic crystalline ice (I c ), 
respectively. Most of the observed 3.07 pm absorption bands in 
dense molecular clouds are best represented by a low-T(<30 
K) deposit (Hagen et al. 1981) and, therefore, by the high- 
density form of amorphous water ice. The frost on interstellar 
grains constitutes the largest reservoir of ice in the unverse — in 
our Galaxy alone, an estimated 10 37 kg — and we therefore 
conclude that most ice in the universe is in this high-density 
amorphous form. 

In the near future, the Infrared Space Observatory will offer a 
unique means to measure the structural state of water ice in 
different astrophysical environments, either directly, from the 
shape of the infrared absorption and emission bands of the 
vibration and bending modes of water, or indirectly, from 
the behavior of trapped impurities in the water matrix. 
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